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ABSTRACT 

We study the possibility that the universe has compact topologies T 3 



T 2 



or S 1 x R 2 using 
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the seven-year WMAP data. The maximum likelihood 95% confidence intervals for the size L of the 
compact direction are 1.7 < L/Lq < 2.1, 1.8 < L/Lq < 2.0, 1.2 < L/Lq < 2.1 for the three cases, 
respectively, where Lq = 14.4 Gpc is the distance to the last scattering surface. An infinite universe is 
compatible with the data at 4.3 er. We find using a Bayesian analysis that the most probable universe 



has topology T 
Subject headings: 



with L/L Q = 1.9. 



1. INTRODUCTION 

General relativity is a local theory and does not pre- 
dict or constrain the global properties of the spacetime 
manifold describing our universe, which have to be con- 
strained through observations. Considerations of non- 
trivial spatial t opology go back to a s early as 1900 by 
Schw arzschild (|Schwartzschildl |l900) (English transla- 
tion (jStewart et al.lll998D b In 1917 de Sitter noted that 
Einstein's equations of general relativity did not con- 
strain the global structure of the spacetime (|de Sitter! 
1917), while in 1924, Friedmann realized that non- 
posit ive curvature spa ces could also have compact topol- 
ogy (jFriedmannl 1 1 9241 ) . A number of other papers were 
published in the follow ing years (for reviews on cos- 
mic topology s e e, e.g. |Lachieze-Rev fc Lumined 119951: 



Luminct 2008; Starkman 1998; Rcboucas & Gomero 



2004 iRoukeTnal [2*000; Blanl oeil fc Roukem a 2000; Levin 



2002[ iBarrow fe Kodamal 120011) , and references therein) 
but the subject remained purely speculative un- 
til very recent accurate measurements of the CMB 
(cosmic microwave backgr ound) anisotr opics by the 
COB E dBennett et al |[l996h and WMAP (I Jarosik et al.l 
l201lHLarson et al.ll2011l:lKomatsu et al.ll2011[ ) satellites, 
and other observational advances in astrophysics and cos- 
mology. 

Various models for the topology of the universe (for 
the classification of different possible topologies see, 
e.g. (|Barrow fc Kodamal I2001D ) have been extensively 
studied recently and compared to the experimental data. 
The two most important ones are the Poincare dodeca- 
hedral space and the 3-torus T 3 ; these models are in best 
agreement with the experimental data. The Poincare do- 
decahedral space arises by slicing the 3-sphere S 3 and 
thus has positive curvature, while the 3-torus is ob- 
tained by slicing infinite Euclidean space R 3 and there- 
fore is flat. Theoretical arguments about quantum cre- 
ation of the universe favor the flat case. Bas ed on 
the W heeler-De Witt equation, Linde has argued (jLindd 
I2004D (see also (|Zeldovich fc Starobinskvl I1984D ) that 
compact flat universes are much easier to create than 
other models, and can naturally provide initial conditions 
for the onset of inflation. Furthermore, Roukema con- 
structed a measure on the set of compact manifolds and 
showed that non-flat models alm ost never occur while 
flat models occur almost certainly (Rouk ema fc Blanloeill 



EoT3) ( see also (|Roukemal [2010) for a discussion of 
the Poincare dodecahedral space versus the 3-torus). 
Following these arguments, we will focus our atten- 
tion entirely on the flat case in this work. Details 
on the Poincare dodecah edral space and experimental 
data analysis are given in (iLuminet et al.|[2003tlRoukemal 
2009t iRoukema et al.l 12004 120081: iCaillerie et al.l 120071: 



Lew fc Roukemal 120 08; Aur ich et al.l 2005b), and other 
topologies in ([Aurich fc Lustid 120101; INiarchou fc Jaffd 
20071: lAurich et all 120061 I2005at ICresswell et al.l I 200C 
Aurich et al.l I2005ct iRoukemal 120051 : lAurich et all 12004 



Riazuelo et al.l 12004 iGomero et al.l l2002al; iBond et all 
2000; iMota et al.ll2005t lGundermannll2005D . 

We analyze three different flat topologies in this work, 
.J?o = T 3 , Jt x = T 2 x R 1 , and .J( 2 = S 1 x R 2 , where 
the subscript denotes the number of non-compact direc- 
tions. We will generically refer to all three cases as a 
torus. Usual flat space is j^^o = K 3 . The topology 
of the 3-torus T 3 can be obtained by identifying the 
opposite edges of a parallelepiped. We only consider 
the simplest case of a rectangular parallelepiped with 
equal side lengths L. This has the highest number of 
symmetries which helps reduce t he computational time. 
Moreover, it has been argued in (We eks et al.ll2004T ) that 
only in well-proportioned spaces is the quadrupole of the 
CMB temperature-temperature correlation function sup- 
pressed compared to the infinite universe. The surpris- 
ingly low observed quadrupole is one of the motivations 
to invoke a compact topology. For our case of the 3- 
torus, well-proportioned means that all three sides should 
be approximately equal. The topologies of the spaces 



T 2 



and S l 



are obtained by compactifying only 



2 or 1 dim ension respective ly. Again, following the argu- 
ment of (jWeeks et al.ll2004T ) and for the sake of simplicity, 
we consider only the case where the compactified dimen- 
sions of T 2 x R 1 have the same size. The size of the 
compact directions will be denoted by L. As L — > oo all 
three manifolds reduce to infinite flat space R 3 . 

Different approaches have been proposed for extract- 
ing information about the topology of the universe from 
the experimental data, the two most important ones be- 
ing the circles-in-the-sky test and the analysis of the 
CMB power spectrum. The basic idea of the circles-in- 
the sky test is that if the global structure of the space 
is smaller than the distance to the LSS (last scattering 



2 



surface), then the LSS will self- intersect in circles, pro- 
ducing correlations between circles with different centers. 
The detection of such circles can reveal the global prop- 
erties of s pace (for deta iled de scription of the method 
see, e.g. (iCornish et all Il99l iMota et all I20T0I [200l 
Levin 2004a b)). The main disadvantage of the method 
is that it cannot be used if the size of the universe is 
bigger than the observable part of it (the distance to 
LSS). The one- year WMAP data has been analyzed with 
this method for signatur es of non-trivial spatial topol- 
ogy (|Cornish et al.l 120041 ) . ruling out the possibility of 
compact spaces with a length-scale smaller than 24 Gpc. 
This limi t has been extended by about 10% by the au- 
thors of (Sha piro Key et aLll2007l) who have also ruled 
out the po ssibility of Poincare dodeca hedral space. The 
authors of (jBielewicz fc Bandavll2011l ) have analyzed the 
most recent seven-year WMAP data with this method, 
putting a lower bound of about 27.9 Gpc on the size of 
the fundamental domain for a flat universe. 

The low-! (i.e. large scale) portion of CMB correlations 
is sensitive to the topology of space, which gives rise to 
another method for detecting the topology. The torus 
preserves the homogeneity of infinite space but breaks 
rotational invariance. This implies that the power spec- 
trum of CMB temperature-temperature correlations does 
not contain all the possible information since the off- 
diagonal elements of the covariance matrix in the spher- 
ical harmonics expansion are non-zero in general, while 
the diagonal elements with equal I and different m values 
are not all equal to e ach other (see sect ion El for more de- 
tails). Moreover, in (|Phillips fc Kogutll2006l) it has been 
argued that the off-diagonal elements contain more infor- 
mation than the diagonal ones if the side length of the 
torus is less than twice the distance to the last scattering 
surface. Therefore, to gain all the possible information 
from the correlations of CMB anisotropies, one has to 
consider the full covariance matrix rather than just the 
power spectrum. 

isly 

used to analyze COBE (Ide Oliviera Costa fc Smootl 
[19951 Ide Oliveira -Costa et al.l Il996fl. on e-year 



The CMB correlation fun ctions have been previously 

/ie 

IT 

WMAP (jPhillins fc Kogutl 120061: IKunz et all I2006T) 

and 
The 



three-year WMAP (jAurich et al 
lower bound on the side length 



120081) data. 
L of T 3 ob- 



tained from COBE data (jde Oliviera Costa fc SmooD 
119951 ) is L > 4.32ft- 1 Gpc at 95% confidence, and 
L > 5.88ft" 1 Gpc at 68% confidence. For T 2 x R 1 
and 5* 1 x R 2 the lower boun d obtained from 
COBE (Ide Oliveira-Costa et alJll99lih is L > 3.0ft- 1 Gpc 
at 95 % confidence. The authors of (|Phillips fc K ogut 
I2006D have obtained higher bounds for T 3 ; L > 1.2 L 
at 95% confidence and L > 2.1 Lq at 68% confidence. 1 
They have also found that the maximum likelihood 
occurs for L — 2.1 L (29 Gpc). Several models of 
tor i with different side lengths have been considered 
in ()Kunz et al.ll2006D with the conclusions L > 19.3 Gpc 
for T 3 and L > 14.4 G pc for 5 1 x M 2 . The main result 
of (jAurich et al.l l2008h is that the 3-torus with volume 
« 5 x 10 3 Gpc 3 (which corresponds to side length of 
17 Gpc) is well-compatible with the WMAP three-year 

1 We will give lengths in terms of Lq = 14.4 Gpc, the dis- 
tance to the last scattering surface. The Hubble length is Hq = 
2. 998/ho = 4.266 (0.703//i ) Gpc. 



data. The WMAP seven-year data has been analyzed 
for detecting signatures of the so-called half-turn space 
(|Aurich fe Lustigl [2010) (the only difference of the 
half-turn space from the 3-torus is that one of the edges 
is turned by 180° before identifying with the opposite 
edge) where the case of the 3- torus is also considered . 
Out of these wo rks, o nly in (jPhillips fe Kogutl I2006D 
and (|Kunz et al.l I2006D has the full covariance matrix 
been analyzed. F or some ear l ier results on these 
topologies see also (|Starobinskyl 119931 : iSokolovl 119931 : 
iStevens et al.lll993l ) and references therein. 

Other methods for experimental detection of non- 
trivial topologies have been proposed and used to an- 
alyze the experimental data. The so-called spati al cross- 
corre lation function method has been used in ()Aurichl 
2008) to analyze the five-year WMAP data for signa- 
tures of a 3-torus topology with a result L = 3.85 Hq 1 . 
Th ey have also specified the o rientation of the torus. 
In ((Bielewicz fc Riazuelol I2008D multipole vectors have 
been used to analyze the five-year WMAP data with 
the conclusion t hat a torus topology is sligh tly preferred. 
The authors of (|Menzies fc Mathews! 120051 ) have looked 
for evidence for a 3-torus topology by the alignment of 
distant objects. They have put a lower bound on the 
side len gth L > 0.9 £p- For some other approaches 
see also dHipolito-Ricaldi fc Gom cro 2005; Marcck i et al 
2005 1 : IGomerol 120031: lHaiian fc Souradeepl 120031: lOpher 



2004; 



2000: 



2002b) 



Dincen et al. 2005; Gomero et al. 2000; Rcboucas 



Gomero et al.l 120011 : IGomerol 120021 : iGomero et al 



In this work we analyze the most recent seven-year 
WMAP data for signatures of the three flat topologies 
of space mentioned above using the full covariance ma- 
trix of temperature-temperature fluctuations. By us- 
ing the symmetry groups of the spaces we construct 
efficient algorithms for the theoretical computation of 
the covariance matrix and the likelihood function us- 
ing that matrix. These algorithms can be used again as 
soon as the high precision da ta from the Planck satellite 
(jPlanck Collaboration! 120061 ) are released. The compu- 
tation of the covariance matrix is done using a mod ified 
version of the CAMB program (jLewis etaDHQOa), as 
discussed in Sec. G2 and that of y 2 and the like l ihood 
using the available WMAP code dJarosik et al.l 120111; 
lLarson et al.ll20lll ; iKomatsu et al.ll201lD ■ We have used 
only the TT correlations in our analysis. Including TE, 
EE, and BB correlations is straightforward, but would 
quadruple the computer time needed, without much im- 
provement in the results since these other correlations 
have much larger errors. 

There have been speculations in the literature 
about the detection of a special direction in the 
CMB map in which the first few multipoles of 
tempera ture-temperature correlations seem to be 
align ed (|Land fc Magueiioll2005l 120071; iRakic fc Schwarzl 
120071 ). This is referred to as the "axis of evil" and is 
given by b = 60°, I = —100° in galactic coordinates. 
The topologies that we consider are not rotationally 
invariant, in particular T 2 x Mr, and S 1 x 
special direction (the infinite one in T 2 



I 2 have one 
I 1 , and the 



finite one in S 



'), so we analyze the case where 



this special direction coincides with the axis of evil, to 
see if the axis of evil can b e explained by one of these 
topologies. The authors of (|Cresswell et al.ll2006D have 
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analyzed the topology S 1 x R 2 with the conclusion that 
it is not the explanation for the multipole alignment. 

The work is organized as follows. In section [2] we 
present a slight generalization of Linde's argument for 
the quantum creation of compact universes. We describe 
the calculation of the covariance matrix and the likeli- 
hood in sections [3] and 2] respectively. We present our 
numerical results in section O and discuss the goodness 
of fit, and maximum likelihood confidence intervals. We 
have done several checks of our analysis, which are given 
in Sec. El The possibility that our results are generated 
by a random fluctuation are analyzed in section [71 where 
we discuss Monte- Carlo skies. The possibility of spuri- 
ous effects due to a small residual CMB dipole in the 
data is investigated in Sec. [8j We summarize in section 
[9l Unless otherwise stated, everywhere in this work the 
side length of the torus is given in units of the distance 
to the last scattering surface Lq. 




a 

Fig. 1. — Plot of the potential U(a) for a positively curved space 
(/3 > 0). The axes are in arbitrary units. 



2. QUANTUM CREATION OF COMPACT 
UNIVERSES 

Consider the standard Einstein-Hilbert action of grav- 
ity minimally coupled to matter. Here we will be only 
interested in "quantizing" gravity, so for the matter por- 
tion we will just consider energy density V without wor- 
rying about where it comes from. Then the action takes 
the form (h = c = 1, M pl = (8ttG)" 1/2 = 1) 



-g(R-2V) . 



(1) 



Following the standard procedure to derive the Wheeler- 
DeWitt equation we use the ADM form of the spacetime 
metric ([Arnowitt et al.ll2008ft 



ds 2 = -N 2 dt 2 + hijidx 1 + JVMi)(daJ + N'dt) 
The action can be rewritten in the form 
S = J d 4 x , 

with 



-^{EijEV - E 2 ) - 2Nv\ 



where 3 R is the 3-curvature of spatial slices, 



E„ 



1 



(2) 
(3) 
(4) 

(5) 



There are numerous possibilities for the spacetime 
manifold and it may be described by infinitely many pa- 
rameters, so to be able to proceed we consider mani- 
folds with finite homogeneous spatial slices which can be 
characterized by one length scale a(t). In other words, 
we assume that locally the manifold is characterized by 
a Friedmann-Robertson- Walker metric while globally it 
can have any finite topology that is compatible with the 
metric. So by a suitable choice of coordinates we get in 
this case 



N = l, Ni = . 
hij = a 2 (t)k 



(6) 
(?) 



where the tensor kij is constant (it only depends on the 
choice of the manifold but does not depend on any of the 



coordinates). Then 



E^EV -E 2 



-G 



(8) 



Since we assumed a homogeneous spatial submanifold 
characterized by single length scale a, by dimensional 
analysis the volume must be proportional to a 3 and the 
curvature to a~ 2 . Namely, 

(9) 



J d 3 x Vh = aa 3 , 







3 R = , 



(10) 

where a and j3 are dimensionless constants that depend 
only on the choice of the manifold. The Lagrangian then 
takes the form 



L = -(a/3- 6aa 2 - 2a 3 V) . 



(11) 



Now we treat a as the dynamical variable describing 
the geometry. The canonical momentum is then 

dL . . 

Pa = ttt = -6aaa , (12) 
oa 

and the Hamiltonian becomes 

H = p a d - L = — !— (-p 2 - 6a 2 (3a 2 + 12a 2 a 4 U) . 

12aa '( 13 ) 

Finally, we canonically quantize, replacing p a by the 
operator — i(d/da) to get for the Hamiltonian 



H 



1 



12aa V da 2 



Qa 2 (3a 2 + Y2o?a*V 



2„4t 



(14) 



Consider now the quantum creation of the universe 
with zero energy. Then the wavefunction of the universe 
"J (a) satisfies the analog of the Schrodinger equation 
with Hamiltonian Eq. (p^|) . which is called the Wheeler- 
DeWitt equation. In this case it takes the form 



— -6a 2 /3a 2 + 12a 2 a 4 V ) *(o) = 0. (15) 
da 2 



The effective potential energy is 

U{a) = 6a 2 Pa 2 - l2a 2 a 4 V , (16) 
and is shown in Fig. [T] which decreases to — oo for large 



4 



a since the second term dominates. However, for small a 
the first term dominates, so for B > there is a po- 
tential barrier from a = to */ B/2V, i.e. the uni- 
verse has to first undergo tunneling before the expan- 
sion can start. This is the reason why the probability 
of quantum creation of positively curved spaces, which 
have B > 0, is thought to be highly suppressed com- 
pared to flat and negatively curved spaces. The action 
for tun neling th rough the barrier is S = a^W /3, where 
ao = v/ B/(2V) is the size of the created universe. The 
tunneling probability is oc exp(— S), and is greater for 
smaller universes; S — > 0, ao — > &s V — > oo. For a flat 
universe, 8 = 0, and the barrier vanishes. 

3. COVARIANCE MATRIX CALCULATION 

Now we turn to the calculation of correlations between 
CMB temperature anisotropies in the flat topologies T 3 , 
T 2 x R 1 , and S 1 x R 2 . Locally they all look exactly 
like the infinite flat R 3 so Einstein's equations and there- 
fore the Friedmann equations are unchanged from the 
infinite case. The calculation for the infinite case is de- 
scri bed in standard textbooks (for a detailed derivation 
see (jDodelsonl 120031 ) ) . so let us briefly summarize that 
calculation and then focus on the differences between 
the infinite and finite universes. Essentially, one has to 
take the Einstein's equations that describe the interac- 
tions between gravity and all of matter and Boltzmann's 
equations for interactions between various types of mat- 
ter (most importantly, electrons and photons) and solve 
for the distribution of photons today given initial condi- 
tions set by inflation. Since the temperature anisotropies 
in the CMB are about five orders of magnitude smaller 
than the background, the calculation is done using per- 
turbation theory around the homogeneous background 
and keeping only first order terms. Then all of the differ- 
ential equations become linear and can be treated easily 
in Fourier space. This is where there is a key difference: 
in an infinite universe the spectrum of the Fourier modes 
k is continuous, while for compactified dimensions the 
spectrum becomes discrete. For a torus with side lengths 
Li, L 2 , and L3 we have k = (ki, ki, fe), 

2?r , 2tt , 2tt 
k\ = —ni, k 2 = —n 2 , k 3 = — n 3 , (17) 
L\ L 2 L3 

where m, n 2 , 713 are integers (the torus is essentially 
a box with periodic boundary conditions). So all of the 
equations in Fourier space remain unchanged, all we have 
to worry about is integrations over k which have to be 
replaced by sums 



d 3 fc 

(27T) 3 



1 



L1L2L3 



(18) 



over the discrete k values in Eq. (fT?)) . The set of points 
Eq. (TT71) will be referred to as the k grid. 

The three cases studied here can be characterized by 
different values for Li. The three-torus T 3 has L\ = L 2 = 
L3 = L, T 2 x R 1 has L\ = L 2 = L, L3 = 00, and finally 
S 1 x R 2 has L\ = L 2 — 00, L3 = L. All three cases 
can be treated in a unified manner by using the integral 



notation, with the understanding that the integral is to 
be replaced by a summation if the corresponding Li is 
finite. 

The first set of summations over k arises when con- 
structing collision terms in Boltzmann's equations. How- 
ever, we will not worry about these integrals for the fol- 
lowing reason. The Boltzmann's equations are important 
only before the decoupling epoch, which corresponds to 
a redshift of about z ~ 1100. The comoving horizon at 
that time was about 50 times bigger than currently, and 
the current bounds on the size of the torus are of the or- 
der of the size of horizon, so at the epoch of decoupling, 
the size of the torus was at least about 50 times bigger 
than the causally connected part. As we will see later in 
section [5l the sums rapidly converge to the correspond- 
ing integrals when the topology scale is around 3 times 
the radius of horizon, which implies that the effects of 
finiteness can be safely ignored for the epoch of decou- 
pling (and before), and k can be treated as a continuous 
variable. All the equations are solved in Fourier space, 
exactly as for the infinite case. 

There is a summation over k when the final answer 
for the temperature fluctuations has to be converted 
from Fourier space back to real space. So let us pick 
up from that point in the calculation. The tempera- 
ture anisotropies 0(n, x) in direction n at a given point 
x (chosen to be our location Xq) are decomposed into 
spherical harmonics 



9(n,x) 



E 



l 

E 



a/ m (x)lj m (n) 



(19) 



l—l m- 



-1 



where the position space ai m coefficients are given in 
terms of the Fourier space temperature fluctuations 
0(n,k) by 



/ S 6 *"/ d " Y im(*M^). (20) 



az m (x) = 



The observed CMB fluctuations are given by the corre- 
lations between the different a; m 's, 



(a/ m (x )a*, m ,(x )) 



(21) 



The correlations between temperature anisotropies in 
fc-space are related to the initial matter power spectrum 

(6(k, n)6*(k',n')) 

where the ma tter power spectrum is defined by 
(jDodelsonll20"0l 

(<5(k)<T(k')) = (2^) 3 5 3 (k-k')P(fc). (23) 

The ratios Q/S on the right hand side of Eq. (|22|) do 
not depend on the initial conditions since the equations 
are linear. All of the information about initial conditions 
is now absorbed into P(k). From Eq. (|20]l. (j2l]t . and ([22} 

we get 



J 



M, 



ImV 



d 3 fc 

(2tt) 3 



P(k) / dn YC m (n) 



8(fc,k-n) 



Q*(fc,k-n') 
' [ ' 5*(k) 



(24) 
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Expanding 0(fc, k • n) into Legendre polynomials 

6(fc,k-n) =^(- J ) i (2/ + l)F i (fc-n)e i (fc), (25) 



and using the identity 

dn P r (fc-n)y im (n) 



47T 

21 + 1 



Su>Yi m (k) 



(26) 



we finally get 

M; m ;/ m < = (47r) 2 (-i) z / x 



d 3 fc 



rP(fc) 



gKfc) er,(fc) 

S(k) S*(k) 



y^(k)Y,, TO /(k). (27) 



We first review the standard result for an infinite uni- 
verse. In this case, the angular part of the integral over 
k in Eq. ([2"7| can be done analytically giving 



Sll'S mm 'Cl 



with 



Ci = - I dkk 2 P(k) 

IT ' 



e,(fc) 



(28) 



(29) 



The derivation remains the same in the finite case ex- 
cept that the k integral must be replaced by the sum 
Eq. ([15]). so instead of Eq. fl57J) we get 

M lmVm , = (4tt) 2 (-»)'*'' x 

' ^ P(k)W-W-YL(m lm ^) . (30) 



5{k) S*(k) 



Now we have to compute a three-dimensional sum 
Eq. (I3U1) instead of a one-dimensional integral Eq. (|29l) 
which requires much more computational time. Also, we 
have to calculate all matrix elements with different I, m, 
I', ml whereas in the infinite case all 1 / 1' or m ^ ml 
(non-diagonal) elements vanish, while the diagonal ones 
do not depend on m. The reason for this is clear. In the 
infinite case, the problem has full rotational invariance, 
so that angular momentum is conserved. In the cases we 
consider, rotational invariance is broken. Even though 
rotational invariance is broken, there is still a large resid- 
ual discrete symmetry group which can be used to sim- 
plify the problem, and reduce the computational time. 
We will refer to this residual symmetry group as G. For 
the T 2 x R 1 and S 1 x R 2 cases, G is the symmetry group 
of a rectangular parallelepiped with two sides equal, the 
tetragonal group D^h with 16 elements, whereas for T 3 , 
G is the symmetry group of the cube, the octahedral 
group Oh with 24 elements. 

The angular part in the sum in Eq. (|30[) can 
be separated (this h as been suggested earlier in 
(|Phillins k K oeut 200fih ) 



Mi 



Iml'm' — 



(4tt) 2 H)^ 



e,(fc) QUk) 



LiL 2 L 3 VV 5 ( k ) S *( k ) 

X £ Yhn&Yvm'fa , (31) 

|k| = fc 

where the first sum is over all the allowed spheres in the 
k grid while the angular sum is over a fixed sphere and 
depends only on the choice of that sphere. 

We can simplify the computation using the discrete 
symmetry group G of the manifolds ^0,1,2- Consider a 



fixed sphere of radius k. If the point (8, </>) of the sphere 
is on the grid, then so are (6, (j> + tt/2), (0,<f> + tt), and 
(9, (f> + Sir/ 2). The angular sum over these four points is 
proportional to 

e i(m'-m)tj) A _|_ e i(m'-m)f _|_ e i{m'-m)-K _|_ e *(m'-m)3f J 

which is unless ml — m is divisible by 4, in which case 
it becomes 4e^ m _m )^. Consider the points (0,<j>) and 
(—9, </>), which both lie on the sphere. Since 

Y lm (-6,(f>) = {-l) l - m Y lm {6,<t>), 

the sum over those two points is unless I + V — m — ml 
is even, but m + ml is even if ml — m is divisible by 
4, so the extra condition we get is that I' — I has to be 
even (this also follows from parity). The eight points 
(±#, (j) + nir/2), n = 0,1,2,3 lie in the eight different 
octants, so the point (9, <fi) can be chosen to lie in the 
first octant. To summarize, the angular sum is nonzero 
only if Z' — Z is even and ml — m is divisible by 4, in which 
case it is equal to 8 times the sum over one octant. Extra 
care is needed for points on the boundary of the octant 
to avoid double counting. 

Consider the points (6, 4>) and (9, tt/2 — <ft) correspond- 
ing to swapping n\ with rii- Taking into account that 
ml — m is divisible by 4, we get 

e i(m'-m)0 + e i(m'-m)(§ -0) = 3 CQS ^ _ m ) 0) ? 

which implies that the angular sums are real. Further- 
more, is also real for even V — I and P(k) and 
Qi(k)/S(k) are real, so the covariance matrix elements 
Mimi'mi are all real implying Mi m Vm' = Afj'm'Jm- Also, 
since Yi_ m (6, (f>) = (— l) m Y^(^, 0) and ml - m is divis- 
ible by 4, we get Mi m Vm> = Mi t - m ,V -m' ■ 

T 3 has more symmetries which can be used to further 
speed up the calculation for this case. For example, the 
sums in Eq. pip over the spherical harmonics are the 
same for all L. Changing L is a rescalling of the allowed 
momenta by 1/L. Thus the angular sum for |k| = k 
for a T 3 of size L is the same as the angular sum for 
|k| = Afc for T 3 of size L/X. Thus the angular sums can 
be computed once, and then used for all values of L. 

Since the calculation of ®i(k)/5(k) is identical to the 
case of infinite fl at universe, we u se the well-known 
CAMB software dLewis et all 120 001 (based on CMB- 
FAST (|Seliak k Zaldarriagalll996[ )) for that part of the 
calculation. Our code is essentially a modification of 
CAMB. It takes the sides of the torus as extra input 
parameters and outputs not only C; but also the com- 
plete matrix M; m ;/ m /. For large values of /, the discrete 
sums over k approach the continuum result, so we only 
use Eq. (|3T|) for I < 30, and use the continuum result for 
Z > 30. The difference between the discrete and contin- 
uum values for Mi m i> m i is less than 0.5% for I = 30. As 
an example, in Fig. [21 we have plotted the ratio of the 
power spectrum G for M\ = T 2 x R 1 with L/L a = 1.9 
to that for infinite space M 3 , where C; has been defined 
as 



Ci = 



1 



1 



21 + 1 



E 



M, 



I nil in 



(32) 



=-/ 



for Jt\ 



T 2 



I 1 . The I = 2 power is reduced by 



7%, and the ratio of power spectra oscillates and rapidly 
approaches unity. The two differ by 0.1% at I = 30. 
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Fig. 2.— Plot of the ratio of C; for = T 2 X I 1 with L/L 
1.9, to that for infinite space R 3 . 



4. LIKELIHOOD CALCULATION 

The matrix M; m ;/ m ' computed as discussed above is 
compared to the experimental data from the 7-year 
WMAP survey. Since rotational invariance is broken, we 
need to vary the orientation of the torus relative to axes 
fixed in space, to find the best fit. We do this by rotating 
the data relative to the torus in computing the likelihood 
function. We specify the orientation by three Euler an- 
gles (<p,9,tp) in the following way. The axes x,y,z are 
fixed in the coordinate frame of the CMB data, i.e. the 
observed universe, and the x' , y' , z' axes are fixed in the 
torus. Start with the CMB-fixed and torus-fixed axes 
aligned. Rotate the torus counterclockwise around the 
z-axis by angle <p, then around the new x-axis by angle 
9, then around the new z-axis by angle ip to get the final 
torus orientation. The angles <p and 9 give the orienta- 
tion of the z-axis of the torus while the angle tp gives the 
orientation of the torus around its z-axis. We can make 
use of the symmetries of our topologies to speed up the 
calculation since various Euler angles can give equiva- 
lent orientations of the torus. Two sets of Euler angles 
(</>, 9, ip) and (</>', 9',ip') are equivalent if 3g <G G such that 



R(^9^) = R{g)R(ct>',9'^') 



(33) 



where R((p, 9, ?p) is the coordinate transformation rota- 
tion matrix corresponding to (<fi, 9, ip) and R(g) is that 
corresponding to the discrete element g. This defines 
an equivalence relation on the set of all possible Euler 
angles. We take a uniform grid on all possible angles, 
then divide that grid into equivalence classes according 
to Eq. (|33|) and take one representative from each class. 
We have scanned over ~ 4000 inequivalent angles. 

Different orientations of the torus were considered 
in t he previous analysis o f first-year WMAP data for 
T 3 ([Phillips fc Kogutll2006[ ). but they only considered a 
uniform grid on the range < cf>, 9, ip < ir/2. Note that 
this does not cover all possible orientations of T 3 . The 
first two angles describe the orientation of the z-axis and 
by their assumption on the range of <p, 9, tp, the z-axis al- 
ways lies in the first octant. However taking into account 
all the symmetries of the cube there are 6 equivalent axes 
that can play the role of the z-axis, the ±x, ±y, and ±z 
axes, while there are 8 octants. In other words, there are 



possible orientations of the cube for which none of the 6 
axes lies in the first octant. 

After choosing a torus orientation, we calculate the 
likelihood in the real space of orientations on the last 
scattering surf ace. For ' N„ p ixels the likelihood function 
is given by 2 (|Dodelsonll2003] l 



5? = 



' — expf-Vc^A) , (34) 



(27r)^/ 2 (det Cyl 
and the \ 2 function by 

x 2 = A T C-'A 



(35) 



where A.; is the vector of pixels and C%j is the covariance 
matrix that also includes the noise. The indices i,j label 
the different pixels, which are in directions n.;j on the 
sky. 

We now describe how to convert the matrix Mi m ii m i to 
the covariance matrix Cy. The t emperature fluc tuation 
measured in a pixel i is given by (|Dodelsonll2003h 

O i= / dn e(n)fli(n), (36) 

where Bi is the beam pattern at the pixel i and is specific 
to the experiment. Usually the beam patterns have the 
same shape for every pixel and are axially symmetric 
around the center of the pixel, as is the case for WMAP, 
so if we denote the direction to the center of the pixel 
by rij then the beam pattern can be decomposed into 
spherical harmonics 



Bi(n) 



(37) 



Using Eq. (fT9f to decompose O(n) into spherical har- 
monics, we get for the theoretical covariance matrix 



Ci 



(e<e i )= tf w m 'B ! ^y im (n 1 )y;;,(n J ), 



Iml' 



(38) 

In computing Cy, we have to vary the orientation of 
the torus relative to the sky. In implementing the Euler 
angle rotation, one can compute the M; m ;< m / matrix in 
the torus-fixed coordinate system, so that it remains un- 
changed as the Euler angles are varied. The pixel direc- 
tions THi are changed to — > R((p, 9, ip) n^. Equivalently, 
one can work in the CMB-fixed coordinate system, and 
rotate the torus, which gives Mi m ii m > transformed by the 
angular momentum rotation matrices, 



A!, 



Iml' 



— )• 



E 



M lnl , n ,D«>(R)D®(R). (39) 



Note that in the infinite universe case, Eq. (|28|) holds, 
and the result Eq. (|38|) simplifies to 



°« - Lu 21 + 1 



BfCiPii&i-Aj), 



(40) 



independent of the rotation R(<p, 9, ip). 

The computation of the covariance matrix using 
Eq. (|38|) is more involved than the infinite case, Eq. (|40"1) , 
so the likelihood calculations require far more computer 

2 To distinguish from length scale L, the likelihood is denoted 
by _S? everywhere in this work. 
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time than the conventional case. dj must be recalcu- 
lated for each set of Euler angles. There are 458403 in- 
dependent elements in Mi m i> m > for 2 < I < 30 of which 
57840 satisfy the I = /'(mod2),m = m'(mod4) condi- 
tion, and 2482 values for each of the indices i and j. The 
slowest step in the computation is evaluating the sums 
on l,m,l',m' in Eq. ([55} for all values of {i,j}- 

An Euler angle rotation of the sky maps points on the 
sphere to rotated points on the sphere. For infinitesimal 
pixels, this corresponds to a reshuffling of the pixels, i.e. 
if pixel i at rij is mapped by the rotation to n., , then pixel 
i — > pixel j . An exact reshuffling of pixels would greatly 
simplify the computation — instead of recomputing Cy , 
one could simply permute the indices on Cy to get the 
transformed matrix. In particular, det C would remain 
invariant under this transformation. 

The WMAP pixels have bee n chosen using the 
HEALPix grid ([Gorski et all 120051) . The pixels are cho- 
sen to lie along lines of constant lattitude, and they have 
equal solid angles. This implies that the spacing of the 
pixels varies as a function of lattitude. As a result one 
cannot treat rotations of the sky as a pixel reshuffling 
transformation. One can approximate the rotations by 
a pixel transformation by mapping the rotated pixel to 
the one closest to it in the HEALPix grid. The likeli- 
hood computed using this method differs from the exact 
result using Eq. (|38[) . and is not accurate enough for our 
purposes. The above approximate relation between rota- 
tions and pixel permutations does, however, explain why 
det C is approximately independent of the Euler angles. 

For a finite universe, one has to use Eq. (J5SJ with 
the value for Mi m i> m i computed as described in Sec. [3J 
The finiteness of the universe only affects the large-scale 
anisotropics, so the difference between the infinite and 
finite cases goes to zero with increasing I. For that rea- 
son we will look only at low-Z portion of anisotropies, 
I < 30, and use the infinite manifold result Eq. (|40p 
for I > 30. We calculate x 2 an d the likelihood _§f us- 
ing a modificatio n of the likelihood code provided by 
the WMAP team dJarosik et al.l[20lll : lLarson et al.ll2011fc 
iKomatsu et al.ll201ll ) as a function of the new parameters 
(L, <fi, 9, ip). Since we are interested only in low-/ effects 
we use the low-resolution portion of the likelihood code. 
We use the experimental data in the exact same form as 
provided by the WMAP team without any further mod- 
ifications. The temperature map used is the smoothed 
and degraded ILC map with the Kp2 mask applied to 
remove the galactic plane and strong point sources. The 
map originally has 3072 pixels, but only 2482 are left 
after the mask. 

Ideally, one would have to do a fit to the experimen- 
tal data varying the four new parameters (L, <p, 9, ip) in 
addition to all the other cosmological parameters. The 
cosmological parameters affect the whole spectrum of 
anisotropies while only the low-Z part of the spectrum is 
affected by the new parameters, so we fix the other cos- 
mological parameters at their best-fit values as give n by 
the seven-year WMAP data (jKomatsu et all 120111 ) and 
only vary the new parameters. The va lues of the cos- 
molo gical parameters that we use are (jKomatsu et al.l 
|20H 100n b h 2 = 2.227, fl c h 2 = 0.1116, n A = 0.729, 
n s = 0.966, t = 0.085, A| J (0.002Mpc- 1 ) = 2.42 x 10~ 9 . 
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Fig. 3. — Plot of In det C/Cf against the Euler angle %f> for fixed 
b,9 for the .^q topology at L/Lq = 1.8, 



5. RESULTS 

Wc have computed the likelihood and x 2 for the three 
cases, J{q = T 3 , Ji x = T 2 x I 1 and 1 2 = S 1 x I 2 , 
for different values of L/Lq as a function of the Euler 
angles. L/ Lq ranges between a minimum value of 0.5 — 
1.0 depending on the manifold, and a maximum value of 
L/Lq = 2.6, in steps of L/Lq = 0.1. By the time L/Lq = 
2.6, the results are almost identical to the flat-space case 
,/#oo = M. 3 . For the statistical analysis discussed later in 
this section, we have used interpolation to construct a 
smooth function of L. 

The relation between likelihood and x 2 is 

- 2 In = x 2 + In det C/C f + In det(27rC» (41) 

where Cf is a fiducial covariance matrix used by the 
WMAP collaboration. Cf is independent of L and the 
Euler angles, and drops out of all likelihood ratios, x 2 
and — 2 In Jzf differ by In det C/Cf (up to an irrelevant 
constant). The likelihood for ^too, three-dimensional flat 
space, will be denoted by Jzfoo, and is Jzfoo = 3573.4 

As noted earlier, for fixed L/Lq, In det C/Cf varies 
weakly with the Euler angles. In Fig. |3j we have plotted 
the variation of In det C/Cf as a function of ip, for fixed 
values of the <p,9, at L/Lq = 1.8. The overall variation 
of In det C/Cf against ip is less than unity. 

— 2 In Jzf (and hence x 2 ) has a strong variation with Eu- 
ler angles at fixed L/Lq. In Fig.|4j we have shown plots 
of the variation of 2 In Jzfoo — 2 In Jzf with Euler angle ip 
for fixed </>, 9 for the topology. Jzfoo is independent 
of the Euler angles. The solid red curve has been chosen 
to have L/Lq = 1.8, and <j),9 values that maximize the 
likelihood at this value of L/Lq. There is a large varia- 
tion of —2 In Jzf with the remaining Euler angle ip, and 
the global minum of -2 In Jzf is 2 In Jzfoo - 21nJSf = -17.2 
at ip/(2w) ~ 0.05. The strong dependence of —2 In Jzf on 
orientation makes it difficult to find the true global min- 
imum of the —2 In Jzf and x 2 functions. We have done a 
scan over all Euler angles with a spacing of 0.057T, to iden- 
tify valleys, followed by a finer scan to find the mininum. 
By comparing our numerical minimum with the next best 
point, we can estimate the uncertainty in our minimum 
—2 In Jzf and % 2 values at less than 0.5. The dashed blue 
curve in Fig. [4]is also for L/Lq = 1.8, but with cp,9 fixed 
at random values, rather than those for which —2 In Jzf 




_i I I I I I I I I I I I I I I I I I I I I I I I 

0.05 0.1 0.15 0.2 0.25 
V/(2n) 

Fig. 4. — Plot of 21n_Sfco — 21ni? against the Euler angle tp for 
fixed <j>, 9 for the topology. The solid red and dashed blue 
curves are for two different values of <j>, 9 at L/Lo = 1.8, and the 
dotted green curve is for L/Lo = 2.2. The solid red curve is for 
L/Lq = 1.8 with <j>, 9 fixed to be the best fit values, the dashed 
blue curve is for L/Lo = 1.8 with <j>, 9 fixed in a random direction, 
and the dotted green curve is for L/Lo = 2.2 with <j>, 9 fixed to be 
the best fit values. 



vs. ip passes through the global minimum. There is still 
considerable dependence as one varies the third angle ip, 
but the dependence is much weaker than for the solid red 
curve. The dependence of —2 In Jz? drops rapidly with in- 
creasing L/Lq. For L/Lo = 2.2, the dotted green curve 
in the figure, the overall variation is about 6.5. 

The plot of x 2 and ~ 21nJz? against L/Lq is given in 
Figs. [5J El [3 for the three cases, J£q = T 3 , J( x = T 2 x R 1 
and ^#2 = -S* 1 x R 2 , respectively. We have plotted the 
maximum and minimum of \ 2 over all possible orienta- 
tions of the torus at each value of L. There is a sig- 
nificant variation in x 2 as a function of orientation, as 
noted earlier, and some orientations are strongly pre- 
ferred over others. In each plot, \ 2 ranges between 
the uppermost and lowermost solid black curves, as one 
varies the orientation of the manifold by varying the Eu- 
ler angles (<fi, 9, ip). For the smallest values of L/Lq, A\ 2 
between the worst and best orientations is 377, 191, and 
156 for ^o,i,2) respectively. As L/Lq increases, the ef- 
fect of a compactified direction decreases. By the time 
L/Lq = 2.6, the fit results are very close to the case of 
the infinite manifold R 3 , and Ax 2 < 4 for the different 
orientations. 

We have been unable to find any pattern to the best- 
fit orientation <fi, 9, tp of the torus as a function of L/Lq. 
We have examined the possibility that the manifolds we 
consider are aligned along the axis of evil. To do this, we 
have chosen the preferred axis of the manifold (the z-axis 
for T 3 , the R direction for T 2 x R 1 and the S 1 direction 
for S 1 x R 2 ) to point along the axis-of-evil direction b = 
60°, I = —100° in galactic coordinates, and allowed for 
arbitrary rotations of the manifold around this direction. 
All the angles are varied with step 7r/100 = 1.8°. This 
gives a subset of all the orientations we have considered, 
and the x 2 range has been plotted as the dashed colored 
curves in the figure. The colored curves lie between the 
black curves (as they must), but they do not lie towards 
the best-fit x 2 ime - This shows that there is nothing 
in our computation that picks out the axis-of-evil as a 
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Fig. 5.— Plot of x 2 and -21njS? against L/L for .£ = T 3 
for different Euler angles (<f>,9,ip). The lower solid black curve 
(solid triangles) is the minimum x 2 (or — -21njSf), and the upper 
solid black curve (open squares) is the maximum x 2 (or — 21n_Sf). 
The lower dashed colored curve (solid triangles) and upper dashed 
colored curve (open squares) are the minimum and maximum x 2 
(or — 21n_S?) with a symmetry axis of the manifold restricted to 
point along the axis of evil. 



C.L. . #o J(\ .#2 
68% 1.471.24 0.92 
90% 1.33 1.02 0.61 
95% 1.24 0.96 0.56 



TABLE 1 

Limits on L/Lo using the x 2 goodness-of-fit test. Values 
less than those in the table are excluded at the 
confidence level given in the first column. 



preferred direction. 

We use a goodness-of-fit test to see whether we can 
rule out the hypothesis that the universe has topol- 
ogy of size L. The minimum x 2 values are x 2 — 
2469, 2467, 2472 at L/Lq = 2.1, 2.1, 2.2 for _# ,i,2, 
respectively. Previous studies have als o found an indica- 
tion o f a dip in x 2 around L/Lq ^2.1 ( Philli ps" fc KogutJ 
2006). Using the computed values of x 2 we have the 
limits given in Table [T] 

We estimate confidenc e inte r vals f or L/Lq us- 
ing likelihood ratio s dCashl 119791 lEadie et"aTI 
l!971f)(|Nakamura et all 120101 §33). Maximum likeli- 
hood confidence intervals are exact if the distribution 
is Gaussian. For non-Gaussian distributions, they have 
an error of order 1/y/n, where n is the number of data 
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Fig. 6.— Plot of x 2 and -2 In Jf against L/L for = T 2 x I 
for different Euler angles (<j>, 8,ip). See the Fig, l5l caption. 



points (|Wiia [19881) . The maximum likelihood value 
is at L/Lq = 1.8,1.9,1.9 for ^#o,i,2, respectively. The 
confidence intervals for L/Lq are determined by using 
AlnJzf, the difference of ln^f f rom the value that max- 
imize s the likelihood function (|Cashl 119791 lEadie et al.l 
1971). Plots of the L/Lq confidence intervals as a func- 
tion of 1 — a, where a is the confidence level, are shown 
in Figs. El QUI The data shows a preference for a finite 
universe with size L/Lq ~ 1.9 corresponding to L ~ 27 
Gpc. The allowed L range extends to L/Lq > 2.6 at 
a confidence level a = 10~ 4 for ^q, 2 x 10~ 5 for 
and 4 x 10~ 3 for Thus the data show evidence for 
a finite universe at a confidence level a — 2 x 10~ 5 for 
the T 2 x M 1 topology. The 95% confidence intervals 
are L/L G [1.7, 2.1] , [1.8, 2.0] , [1.2, 2.1] for .-#04,2, 
respectively. 

We have scanned over ~ 4000 different orientations for 
each value of L, with a finer scan near the minima, so 
that the error in x 2 and —2 In ££ is < 0.5. The difference 
in — 21nj£? between its value at L — ¥ 00 and its minimum 
value (which occurs at L/Lq = 1.9 for is 20.4, which 
is well outside possible numerical errors. Note that the 
main numerical uncertainty is finding the true minimum 
of — 2 lnJ?f for finite values of L. The minimum value of 
— 21n«if has been determined with an accuracy < 0.5. 
The actual difference in likelihoods between finite and 
infinite L can only be greater than what we have found. 
There is an indication that a finite universe fits the data 
better than an infinite one. However, the "standard" 5(7- 
criterion for a discovery, corresponding to a confidence 
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Fig. 7.— Plot of x 2 and -21njSf against L/L for Jti = S 1 X R 2 
for different Euler angles (<j>,6,ip). The two colored curves are 
almost on top of each other. See the Fig. [5] caption. 




Fig. 8. — Plot of the L/Lq confidence interval as a function of 
the confidence level for T . 



level a = 5.7xl0 7 , includes the value L = 00. 

The Euler angles for the best fit case ^#1 with L/Lq = 
1.9 are (<f> = 21° ±2°, 9 = 53° ± 2°,V = 61° ± 2°) 
which corresponds, for the infinite direction, to (b — 
37° ± 2°, I = 291° ± 2°) in galactic coordinates and (a = 
182°±2°, 8 = -25°±2°) in J2000 equatorial coordinates. 
This is close to the direction (b = 30°±2°,Z = 276° ±3°) 
of the velocity of the Local Group inferred from the CMB 
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Fig. 9. — Plot of the L/Lo confidence interval as a function of 
the confidence level for T 2 X K 1 . 




Fig. 10. — Plot of the L/Lq confidence interval as a function of 
the confidence level for S 1 X R 2 . 



dipole (jLineweave JT991 . We discuss the possibility that 
our signal is due to a dipole contamination in Sec. [HJ 

We can compare the probabilities that the universe 
has topology ^#o,i,2,3- We use the Bayesian prior that 
there are four discrete choices, with L/Lq = 1.8, J%\ 
with L/Lq = 1.9, ^#2 with L/Lq = 1.9, or the infinite 
case ^#oo, each of which is equally probably. Since we 
are comparing four discrete cases, there is no ambiguity 
due to choice of measure in choosing an equiprobable 
prior. Then using the likelihood ratios gives the posterior 
probabilities 



o) = 0.04, 
= 0.96, 
2 ) = 2 x 1Q- 



p = 1.3 x 10~ J . 

The probability of the infinite universe j t 
and the most probable topology is JK\ 
two dimensions compactified, and one infinite. 



is very small, 



with 



5.1. Fisher Information 

The Fisher information can be used to compute the 
variance V of the length L determined using the maxi- 
mum likelihood method. The Fisher information is given 
by 



V~ 




TrC 



dC- 



dL 



l c BC 



dL 



ir ° dL C dL 



(42) 



Using the covariance matrix C for the toplogy with 
L/Lq = 1.9 gives 



y- 1 = 3.3 x 10 3 



(43) 



so that the error estimate for L/Lq is \/V = 0.017. The 
Fisher information error Eq. (|43p corresponds to using a 
quadratic approximation to the likelihood function about 
its minimum to determine the error, and gives a smaller 
error than that obtained earlier using the exact likelihood 
function. 

6. CHECKS 

We have been unable to find a simple explanation for 
the better fit due to a finite topology. However, there are 
some possibilities which we can test. 

The measured cosmic microwave background 
anisotropy has a smaller value for the quadrupole 
power C2 than the theoretical expectation value. There 
is a large cosmic variance in C2 , so this is not a discrep- 
ancy between theory and experiment. Fig. [2] shows that 
the predicted value of C2 for a finite universe is reduced 
from the infinite universe value. The greater likelihood 
for a finite universe is not due to lowering the value of 
C%. We have checked this by determining the likelihood 
using M; m ;/ m / for the finite case, but with the I = V = 2 
values replaced by their values for the infinite universe. 
For jtfx w hh L/Lq = 1.9, -21nJz? increases by 0.43, 
which is much less than the 20.4 difference in — 21nJ§? 
from the infinite universe. 

As another test, we have computed the likelihood for 
the Jt\ topology for L/Lq = 1.9 by using a truncated 
Mi m i/ m > matrix. The truncated matrix is constructed by 
using Mi m ii m i for the finite topology for 5 < I, V < 20, 
and using M/ m // m / for the infinite universe, i.e. M; m // TO / = 
CiSu'S mm i, for I and or V outside this range. A plot of 
the likelihood as a function of the Euler angle -0 for this 
truncated Mi m i' m i is plotted as the dashed blue curve in 
Fig. QT] This can be compared with the likelihood curve 
using the full Mi m i' m i for the finite topology, shown as 
the solid red curve. The dip in the likelihood difference 
to —20.4 is the signal that the finite topology is a better 
fit than the infinite universe. The plot for the truncated 
matrix is similar to that for the full matrix, except that 
the small- angle fluctuations have been smoothed out, as 
is to be expected since higher I terms have been dropped. 
Note that the dip in 2 In J?^ — 2 In ££ is very similar in 
both cases, and the minimum of 2 InJzfoo— 2 In Jz? is nearly 
the same. This shows that the effect we find is not due 
to the low-Z modes (quadrupole, octupole), and is also 
not an edge effect as a result of only using I < 30 in 
the compuation. For 5 < I, I' < 20, Mi m i> m > has 86736 
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Fig. 11.— Plot of 21n.S? o - 21ni? against the Euler angle if} 
for fixed <f>, 6 for the topology for L/Lo = 1.9. The solid red 
curve uses the full matrix M; mi / m / , the dashed blue curve uses the 
matrix truncated to 5 < I, I' < 20, and the dotted green curve uses 
the matrix M im ;/ m /<5;;/, retaining only the part diagonal in I. 

elements of which 11056 satisfy the / = V (mod2),m = 
m! (mod 4) condition and are non-zero. 

The off-diagonal elements in Mi m i' m ' are important for 
the calculations. We have also computed the likelihood 
by retaining only the elements which are diagonal in I, i.e. 
using Mimi'm'Su'. This drops the elements in Mi m Vm'&XV 
which are off-diagonal in I, while retaining the elements 
which are off-diagonal in m for a given I. The likelihood 
with this matrix is the dashed green curve in Fig. QTJ 
With this matrix, the likelihood deviates much less from 
the infinite universe, and the dip near i/j/(2ir) 0.17 is 
much less pronounced. 

7. MONTE-CARLO SKIES 

The results of the previous section were obtained us- 
ing a likelihood analysis of the WMAP7 data. One can 
study whether the better fit of a finite topology is due 
to a statistical fluctuation. Since the big-bang is not a 
repeatable experiment, this must be done by generating 
random Monte-Carlo data for the pixels Aj, and redo- 
ing the analysis for this Monte-Carlo data. To actually 
do this numerically is beyond the computing power we 
have available. Luckily, for the problem at hand, we can 
analyze the Monte- Carlo problem analytically. 

Assume that the pixels Aj are generated by the co- 
variance matrix Coo for an infinite universe, so that the 
probability distribution is 



P (A): 



1 



Vdet(2^Coo) 



exp 



1 



A T C~ 1 A 



(44) 



The likelihood function computed using A, and covari- 
ance matrix C (of a finite universe) is 



21ogJ^ = A T C _1 A + In det (277(7) , 



(45) 



and the likelihood constructed using the covariance ma- 
trix C^o of the infinite universe is 



Let 



21og^ 00 -A T C 00 1 A + lndct(27rC 00 ). (46) 



ft = (-2 log J?,*,) - (-2 log if) , (47) 



be the difference of the two log-likelihoods. In our anal- 
ysis, we found h = 20.4 > 0, so that the finite universe 
was more likely than the infinite universe. The average 
value of h over Monte-Carlo data can be computed using 
Eq. flH) and Eq. ([iT]). The two-point function is 

(A i A j ) = (C 00 ) ij , (48) 

so that 

(h)=N — DrC^Coo + lndet(C* 00 ) - lndet(C) .(49) 

It is convenient to define the symmetric matrix 

S = C^C-'CU 2 , (50) 

which is a positive matrix since C and Coo are positive 
matrices, and has eigenvalues Sj > 0. In terms of S, 



(ft) =N- TrS + lndetS 
= ^[l-Si + lnsi] . 



(51) 



The function 1 - 
s = 1. Thus 



s + ln s < with its maximum at when 



(h) < 



(52) 



and (ft) = only if S = 1, i.e. C = Coo- This gives the 
intuitively obvious result that the best fit for data gener- 
ated with covariance matrix C M is, on average, given by 
fitting using the same covariance matrix Coo ■ Any other 
covariance matrix C used for fitting, on average, gives a 
lower likelihood. 
If instead of Eq. (|47|) we had used the difference of x 2 , 



h = 2 — 2 



then 

(h x )=N -TrC^Co 



N -TrS 



(53) 



(54) 



and (h x ) could have either sign, since Sj > 0, but need 
not be smaller than 1. For example, a simple rescaling, 
C = XCoo, with A — > oo can always make x 2 ~ > 0) its 
minimum poissible value. This option is eliminated for 
likelihood because of the det(27rC) term. 

Using for C the best-fit topology -#i with L/L = 1.9, 
we find the numerical values 



N = 2482 , 
In det(C 00 ) = 1097.8 + In detC/ , 
In dot (C) = 1082.73 + lndet C f 
Tr C- 1 C ao = TrS = 2516.6, 



so that 



(ft) = -19.5. 



(55) 



(56) 



This differs from the value we find of h = +20.9 by Ah = 
h — (h) = 40.4. The probability that Ah is a statistical 
fluctuation can be determined by computing the variance 
of h using the four-point function 

(AjAj-AfcAj) = (q„)y (Coo) kl + {Coo) lk {Coo)^ 

+ (Coo) u (Coo) jk , (57) 
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to obtain 

((Aft) 2 ^} = 2N + 2TrC- 1 C 00 C- 1 C 00 - 4 (Tr C^C*,) 
= 2Tr (l-C- 1 C 00 ) 2 = 2Tr(l-S') 2 . (58) 



In our case, 



so that 



Tr C~ 1 Coo C~ 1 Coo = 2610.0 , 



(Aft) 2 ) = 117.6= (10. 



(59) 



(60) 



Our observed value of Aft = 40.4 is 3.7 tr away from 
the mean, so the probability that a fluctuation gives h 
larger than or equal to our observed value is 1.1 x 10~ 4 , 
assuming a normal distribution. 

While the distribution of the data Aj is Gaussian, the 
distribution of the likelihood difference h is no longer 
Gaussian. We can also compute higher order connected 
correlation functions of ft, 

((Aft)'') c = 2 r (r-1)! Tr (l - C^C^Y 

= 2 r (r-1)! Tr {1 - S) r , (61) 

from the generating function 

log(e AA,l ) = -ATr(l-C- 1 Coo) 

-^Trln [l-2A(l-C" 1 C 0O )] 

= -ATr (1 - S) - ^Tr In [1 - 2A(1 - S)} , 



so that 



(62) 



((Aft) 3 ) = 8Tr (l-C-^oo) 3 
((Aft) 4 ) c = 48Tr (l-C^^oo) 4 , (63) 



where the fourth-order correlation is 

\4\ //a i,\4\ 



<(Aft) 4 ) = ((Aft) 4 ) c + 3<(Aft) 2 y 



(64) 



(Aftf\= -489.0, 
(Afr) 4 \ =4303.7. 



in terms of the connected correlation. The mean value 
(ft), and all the connected correlation functions ((Ah) r ) c 
are of order N, the number of data points. Thus the 
relative correlation ((Aft) r ) / (h) r is of order iV 1 ~ r . 
The numerical values for our case are 



(65) 



We can get a better estimate of the probability that 
h = 20.4 is due to a statistical fluctuation by using these 
higher order moments. We have fit (ft) and ((Ah) r ), 
r = 2, 3, 4 to a probability distribution 

p(h) = p esq) [-(ft - ft ) 2 - c 2 {h - h Q ) 2 

~c 3 (h~h a ) 3 -c A (h-h )' i ], (66) 

and found using this distribution that the probability 
that ft — (ft) > 40.4 is 10 -6 , which is smaller than the 
value obtained earlier using a normal distribution for ft. 



8. DIPOLE CONTAMINATION 

The symmetry axis of points in the direction 
(b = 37° ±2°,l = 291° ± 2°), which is close to the direc- 
tion of the v elocity of the Loca l Group (b = 30° ±2°, I = 
276° ± 3°) (|Lineweaverlll996D . The CMB has a large 
dipole asymmetry of 3.358±0.001±0.023 mK in the direc- 
tion (b = 48.05° ± 0.11°, I = 264.31° ± 0.2°) (jLineweaverl 
Il996[ ) . Suppose that the data is contaminated by a dipole 
contribution that has not been properly subtracted out. 3 
Could a residual dipole explain the results we have 
found? 

To study the effect of a residual dipole, assume that 
the observed pixels are 

A° bs = Aj +p • &i = Ai + <k, di = p ■ n 4 (67) 

where Aj are the true fluctuations given by the distribu- 
tion Eq. (|44[) and p is the residual dipole contamination 
in the data. Then Eq. (|45|46|) are replaced by 

- 2 log = (A + df C' 1 (A + d) + In det(27rC) , 
-2 log = (A + df C^ 1 (A + d)+ lndct(27rC) . 

(68) 



From these, we find 

(ft)=jV — TrC~ 1 C 00 + lndet(Coo 
+ (d T C^d) - (d T C~ 1 d) , 
(Aft 2 ) = 2A^ - M T C- 1 d + 4d T C 00 1 d 



- lndet(C) 

4d T C~ 1 Co C- 1 d 
(69) 



-2Hr(7- 1 C QO C- 1 C QO - 4 (TyC^C^ 



Dipole contamination produces a systematic shift in h 
from its value in Eq. (j49]) given by the (d T C^ Q 1 d) — 
(d T C~ 1 d) terms, which can be written as 



(d T C- l d) - (d T C^d) = PaP0 D a0 



(70) 



in term of the components p a — {p x ,Py,Pz) of the dipole. 
We find 



D = 



-1.22 

-0.004- 

-0.114 



-0.004 -0.114 
-0.0281 -0.003 
-0.003 -0.0127 



mK 



(71) 



The largest eigenvalue is — 1.23mK -2 . To get a shift in 
— 21nJzf of 20.4 requires a dipole contamination |p| of 
around 4mK -2 . This is larger than the observed dipole, 
and several h undred times the quoted uncertainty in the 
CMB dipole JI ineweave"rlll996h . and is excluded. 

9. CONCLUSIONS 

We have analyzed the possibility that the universe 
has compact topologies = T 3 , = T 2 x I 1 and 
■Mi = S 1 x R 2 using modifications of the available C AMB 
and WMAP 7-year likelihood codes. The maximum like- 
lihood 95% confidence intervals are 1.7 < L/Lq < 2.1, 
1.8 < L/Lq < 2.0, 1.2 < L/L < 2.1 for the three cases, 
respectively. Using the Bayesian analysis discussed ear- 
lier, we find that the most probable universe has the 
compact topology . An infinite universe is compati- 
ble with the data at a confidence level of 2 x 10~ 5 (i.e. 
4.3 a). We find no evidence of a preference for the axis- 
of-evil direction. The improved fit for a finite universe 

3 This possibility was suggested to us by B. Keating. 
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is not due to the lowered prediction for the quadrupole 
anisotropy; this accounts for only a small fraction of the 
increase in likelihood. 

It would be useful to investigate whether any system- 
atic effects in the WMAP data introduce effects with 
cubic symmetry that can mimic the effects of a torus 
topology. The best fit results do not pick out any special 
orientation for the torus, such as a torus with symme- 
try axis perpendicular to the galactic plane, that might 
lead to systematic effects that lead to a fake signal. Pix- 
elization of the data using the HEALPix grid also should 
not introduce cubic symmetry terms along an axis not 
aligned with the galactic pole. The best fit results for 



^#1 have symmetry axis which is near (<- 
tion of the Local Group velocity. 



10°) the direc- 
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